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Scaling behavior of the order parameter and its conjugated field in an absorbing phase 

transition around the upper critical dimension 
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We analyse numerically the critical behavior of an absorbing phase transition in a conserved 
lattice gas in an external field. The external field is realized as a spontaneous creation of active 
particles which drives the system away from criticality. Nevertheless, the order parameter obeys 
certain scaling laws for sufficiently small external fields. These scaling laws are investigated and 
the corresponding exponents are determined in various dimensions (D — 2, 3, 4, 5). At the so-called 
upper critical dimension D c = 4 one has to modify the usual scaling laws by logarithmic corrections. 

PACS numbers: 05.70.Ln, 05.50.+q, 05.65.+b 



I. INTRODUCTION 



Recently Rossi et at. introduced a conserved lattice gas 
(CLG) with a stochastic short range interaction that ex- 
hibits a continuous phase transition from an active state 
to an absorbing non-active state at a critical value of 
the particle density ffl. The CLG model is expected 
to belong to a new universality class of absorbing phase 
transitions characterized by a conserved field. Similar to 
the well known universality hypothesis of directed per- 
colation |, [| the authors conjectured that "all stochas- 
tic models with an infinite number of absorbing states 
in which the order parameter evolution is coupled to a 
nondiffusive conserved field define a unique universality 
class" [jy. In order to check this hypothesis the scal- 
ing behavior, i.e. the values of the corresponding critical 
exponents, has to be determined. The order parameter 
exponent as well as the exponent of the order parame- 
ter fluctuations were determined in Q. Furthermore a 
modified CLG model with random particle hopping was 
introduced in || which mimics the mean-field scaling be- 
havior of the system. 



In this work we investigated for the first time the CLG 
model in an external field. The external field is con- 
jugated to the order parameter, i.e. it is realized as a 
spontaneous creation of active particles. Of course a 
spontaneous creation of active particles destroys the ab- 
sorbing state and thus the absorbing phase transition at 
all. But analogous to ferromagnetic phase transitions 
the order parameter can be described as a generalized 
homogeneous function of the control parameter and of 
the conjugated field in the critical regime. This scaling 
behavior is examined below (D = 2,3), above (D = 5), 
and at the so-called upper critical dimension (D = 4). 



II. MODEL AND SCALING BEHAVIOR 

We consider the CLG model on D-dimcnsional cubic 
lattices of linear size L. Initially one distributes randomly 
N = pL particles on the system where p denotes the par- 
ticle density. In order to mimic a repulsive interaction 
a given particle is considered as active if at least one of 
its 2D neighboring sites on the cubic lattice is occupied 
by another particle. If all neighboring sites are empty 
the particle remains inactive. Active particles are moved 
in the next update step to one of their empty nearest 
neighbor sites, selected at random. In the steady state, 
the system is characterized by the density of active sites 
p a . The density p a is the order parameter of the ab- 
sorbing phase transition, i.e., it vanishes if the control 
parameter p is lower than the critical value p c . In the 
thermodynamic limit the order parameter scales for zero 
field and for p > p c as 



p*(6p,h = 0) ~ Sff 



(1) 
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with Sp = p/ p c — 1. 

In this work we consider the CLG model in an exter- 
nal field h conjugated to the order parameter, i.e., the 
field spontaneously creates active particles. Clearly the 
particular implementation of the spontaneous creation of 
active particles has to obey the particle conservation of 
the CLG model which is believed to be a relevant pa- 
rameter, determining the universality class. Instead of 
a creation of additional particles we realize the external 
field by particle movements which do not change the to- 
tal number of particles. Therefore we choose randomly 
L D h particles on the lattice. Each of these particles is 
then moved to one of its empty next neighbors. In this 
way inactive particles may be activated and the number 
of active sites is increased. 

In our simulations we start with randomly distributed 
particles. The system is updated according to the above 
rules in the following way. One lists all active sites and 
updates these sites in a randomly chosen sequence. Then 
one moves L D h randomly chosen inactive particles to one 
of its empty next neighbors. In the case that these move- 
ments create active particles, these particles are added to 
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FIG. 1: The density of active sites p a as a function of time 
(number of update steps) for two different values of the ex- 
ternal field h. After a certain relaxation time the density of 
active sites fluctuates around a well defined value. For too 
small fields finite-size effects occur which results in strong 
fluctuation events (see lower curve). 

the list of active particles and will be updated in the fol- 
lowing update step. Thus one update step contains both, 
the update of active sites and the additional movements 
of inactive particles which mimics the external field. 

In Fig.jl] we plot the density of active sites vs. time 
(number of update steps) at p c for two different values 
of the driving field. After a certain relaxation time the 
system reaches a steady state where the density of ac- 
tive sites fluctuates around the average value (p a {5p, h, t)) 
which is interpreted as the order parameter p a (5p, h). De- 
creasing the external field below a certain value finite-size 
effects occur which results in strong fluctuations. Two of 
these strong fluctuations events can be seen in the lower 
curve of Fig.g. The origin of these effects is that the 
number of particle movements (L D h) create too few ac- 
tive particles and the system tends to the absorbing state. 
As a results the measured order parameter is shifted to 
lower values. To avoid these finite-size effects we increase 
the system size before these fluctuation events occur. 

The spontaneous creation of active particles destroys 
the absorbing state and thereby the phase transition it- 
self. Although the external field drives the system away 
from criticality the order parameter obeys certain scaling 
laws for sufficiently small fields. At the critical density 
the order parameter is expected to scale with the field as 
(see for instance Q) 

p a (6p = 0,h) ~ hU a . (2) 

As usual in critical phenomena the order parameter is 
assumed to be a generalized homogeneous function of the 
control parameter dp and the applied field h 

P& {5p,h) = Af^A-^U^) (3) 



with the scaling function f. Choosing SpX^ 1 ^ 13 = 1 at 
zero field one recovers Eq. (Q) whereas h A~ <T /' 3 = 1 leads 
to Eq. (||) at the critical density, respectively. 

Furthermore we consider the fluctuations of the order 
parameter 

A Pa (Sp,h) = L D [( Pa (6p,h 7 t) 2 ) - ( Pa (5p,h,t)) 2 ]. 

For zero field the fluctuations are known to diverge ap- 
proaching the critical point [Q 

Ap a (Sp,h = 0) ~ 6 P ^'. (5) 
The fluctuation exponent 7' fulfill the scaling relation Q] 

i = v ± D - 2(3, (6) 

where the exponent i> ± describes how the spatial correla- 
tion length diverges at the transition point. In the critical 
regime we assume that the fluctuations obey the scaling 
ansatz 

A Pa (Sp,h) = A 7 ' d(5 P X,hX°) (7) 

Setting 5p\ — 1 one recovers Eq. (||) for h = 0. 

Analogous to equilibrium phase transitions the suscep- 
tibility is defined as the derivative of the order parameter 
with respect to the conjugated field 

d 

x{Sp,h) = —p^Spih) 

= X 1 -*'? cVpX-W.hX-*'*). (8) 

Approaching the transition point the susceptibility di- 
verges for zero field as 

X (Sp,h = 0) ~ 5p-\ (9) 

This result can be recovered from Eq. (|J) by setting 
SpX 1 ^ = 1 for h = and one gets the scaling relation 

7 = u - P (10) 

which corresponds to the well known Widom equation of 
equilibrium phase transitions. Using this scaling relation 
we calculate in the following the value of the suscepti- 
bility exponent 7 from the obtained values of p and a. 
Notice that in contrast to the scaling behavior of equi- 
librium phase transitions the non-equilibrium absorbing 
phase transition is characterized by 7 7^ 7'. 

III. BELOW THE CRITICAL DIMENSION 
A. D=2 

At the beginning of our analysis we consider the two- 
dimensional CLG model at the critical density p c for var- 
ious values of the driving field h where the value of the 
critical density is obtained from B . The field dependence 
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FIG. 2: The order parameter p a as a function of the field h 
at the critical density p = p c for D = 2. The dashed line 
corresponds to a power-law behavior according to Eq. (0) with 
an exponent 0.286 ± 0.001. 



of the order parameter p a (Sp — 0,h) is shown in a log- 
log plot Fig. |[ Approaching the critical point (h — > 0) 
the order parameter vanishes algebraically in agreement 
with Eq. (||) . A regression analysis of the data yields the 
estimation (3/a = 0.286 ± 0.001. Using the value of the 
order parameter exponent (3 = 0.637 ± 0.009 B we get 
cr = 2.227 ± 0.032. 

In the following we analyse the order parameter as a 
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FIG. 3: The scaling plot of the order parameter p a for the 
two-dimensional model. The data are rescaled according to 
Eq. (|ll|) with 13 = 0.637 and a = 2.227. The inset displays 
the unsealed data, i.e., the order parameter p a is plotted as 
a function of the density p for different values of the field h. 
The order parameter is driven away from the critical point 
(p c ,0) with increasing field. The dashed line corresponds to 
the zero- field behavior. 



FIG. 4: The scaling plot of the order parameter fluctua- 
tions Ap a for the two-dimensional model. The data are 
rescaled according to Eq. ( |l2| ) with a = 2.227 and 7'/ " = 
0.18. The inset displays the unsealed data, i.e., the fluctua- 
tions are plotted as a function of the density p for different 
values of the field h. Approaching the transition point (h — > 
at p = p c ) the peak of the fluctuations diverges. 

function of the control parameter bp for different fields 
from h — 10~ 5 up to 2 10 -4 . The applied field results in a 
rounding of the zero-field curve, i.e., the order parameter 
increases smoothly with the control parameter for h > 
(see inset of Fig. . According to the scaling ansatz of 
the order parameter [Eq. (j^)] we choose h \~ a ^ = 1 and 
get the scaling form 

Pa (5p,h) = hP/° fidph- 1 '",!). (11) 

Thus plotting p a h~^^ a as a function of Sph" 1 ^ the 
curves for different values of the driving field have to 
collapse onto the scaling function f. Using the above de- 
termined values of (3 and cr one gets an excellent data 
collapse which is shown in Fig.||. 

Next we consider the scaling behavior of the order 
parameter fluctuations Ap a . The fluctuation data are 
shown for different values of the external field in the inset 
of Fig. For finite fields the fluctuations display a peak. 
Approaching the transition point (h 0) this peak be- 
comes a divergence signalling the critical point. In order 
to analyze the scaling behavior of the fluctuations we use 
Eq. (0) and set h A CT — 1 which yields 

A Pa (Sp 7 h) = bri'l" diSph- 1 '*,!). (12) 

Since p c and a are already determined one can estimate 
the exponent 7' by varying this exponent until one ob- 
serves a data collapse of the different fluctuation curves. 
The best result is obtained for j' /cr — 0.18±0.02, leading 
to 7' = 0.402±0.045, and the corresponding plot is shown 
in Fig{| The value agrees with 7' = 0.384 ±0.023 M ob- 
tained from a regression analysis according to Eq. (@) . 



4 



10" 1 : 



L= 32 
•L= 64 
■ L=128 
■ L=256 



10' 



x 



D=3 



10° 



10"' 



10' 



10 10 
h 



10° 



10" 



0.20 



0.15 



0.10 

a. 

< 

0.05 




0.00 <= 
-12.0 



-4.0 



0.0 



4.0 8.0 



5p ft" 



FIG. 5: The order parameter p a as a function of the field h 
at the critical density p = p c for D = 3. The dashed line 
corresponds to a power-law behavior according to Eq. (0) with 
an exponent 0.403 ± 0.004. 



Additionally the fluctuation exponent 7' can be esti- 
mated via the scaling relation Eq. (||) . Using the estima- 
tion v ± = 0.78 ± 0.08 we get 7' = 0.286 ± 0.161 which 
again agrees with our result. 



B. D=3 

Analogous to the two-dimensional case we first deter- 
mine the exponent a for D = 3. The field dependence 
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FIG. 6: The scaling plot of the order parameter p a for the 
three-dimensional model. The data are rescaled according to 
Eq. (0) with 13 = 0.837 and a = 2.075. The inset displays 
the unsealed data, i.e., the order parameter p a is plotted as 
a function of the density p for different values of the field h. 
The dashed line corresponds to the zero-field behavior. 



FIG. 7: The scaling plot of the order parameter fluctua- 
tions Ap a for the three-dimensional model. The data are 
rescaled according to Eq. ( |l^ ) with a = 2.075 and j' /cr = 0.1. 
The inset displays the unsealed data, i.e., the fluctuations are 
plotted as a function of the density p for different values of 
the field h. 



of the order parameter is plotted in Fig.|5|. A regres- 
sion analysis yields f3/a = 0.403 ± 0.004 which leads to 
the estimation a = 2.075 ± 0.043 if one uses the value 
(3 = 0.837 ±0.015 J§. 

In order to investigate the scaling behavior of the order 
parameter we simulated the CLG model for field values 
from h = 2 10~ 5 up to h — 2 10~ 4 . The corresponding 
curves as well as the zero-field order parameter curve is 
shown in the inset of Fig. ^. According to Eq. ( |TT|) we plot 
the rescaled order parameter in the same figure. Using 
/3 = 0.837 and a — 2.075 we get again and excellent data 
collapse of the different curves. 

Furthermore we analyse the scaling behavior of the 
fluctuations Ap a (5p, h) for the three-dimensional CLG 
model. The corresponding data are presented in the in- 
set of Fig. [?]. Similar to the two-dimensional case one 
tries to obtain a data collapse of these curves by vary- 
ing the exponent Y/a [Eq. (Jl2[)] . The best result is ob- 
tained for i jo = 0.10 ± 0.02 (see Fig.g) which leads to 
7' = 0.208 ±0.042. 

Previous investigations of the fluctuations Ap a for 
zero-field showed that Ap a diverges for Sp — > 0. But the 
numerical data could be interpreted either as a power-law 
divergence [Eq. (||)] with an exponent 7' = 0.18 ± 0.06 
or as a logarithmic growth Our new results show 
that the fluctuations of the three-dimensional CLG model 
diverges algebraically at the transition point with the 
exponent 7' = 0.208 ± 0.042. Additionally, the above 
results allows to estimate the correlation length expo- 
nents via the scaling relation Eq. (|^) and one obtains 
V\ = 0.627 ± 0.027. 




FIG. 8: The order parameter p a as a function of the field h at 
the critical density p = p c for D — 4. According to the scal- 
ing ansatz Eq. ( |15[ ) we plot p a /i -1//2 as a function of | In 
The expected asymptotic behavior p a /i _1//2 = const | In /i| E 
is observed for E = 0.41 (dashed line). The inset displays 
the original data. The dashed line corresponds again to the 
ansatz Eq. with E = 0.41. 



FIG. 9: The scaling plot of the order parameter p a for the 
four-dimensional model, i.e., at the upper critical dimension. 
The data are rescaled according to Eqs. ( |l6|Jl7| ) using the 
mean-field values j3 = 1 and a = 2. The inset displays the 
unsealed data, i.e., the order parameter p a is plotted as a 
function of the density p for different values of the field h. 
The dashed line corresponds to the zero-field behavior. 



IV. AT THE CRITICAL DIMENSION 

At the upper critical dimension D = 4 the scaling 
behavior of the CLG model is affected by logarithmic 
corrections Q, i.e., the scaling ansatz Eq. (|J) has to be 
modified. Motivated by the scaling behavior of the Ising 
model we assume that the order parameter obeys in lead- 
ing order the ansatz (see Appendix) 



Pa(6p, h) 



(13) 



A | In \\ l ?(5p\- 1 l \ In A| b , h \^ cr ^\ In A| 5 



where the exponents (3 and a are given by the corre- 
sponding mean-field values /3 = 1 and er = 2, respectively. 
Thus, for zero field the asymptotic scaling behavior of the 
order parameter obeys 



p a (8p,h = 0) ~ SpllnSpf 



(14) 



with B = b + I. This behavior was already observed for 
the CLG model with the logarithmic correction exponent 
B = 0.24 §. 

According to the above scaling ansatz the asymptotic 
field dependence of the order parameter at the critical 
density is given by 



p a {5p = O.h) ~ h 1 ' 2 |ln/i| 5 



(15) 



with £ = s/2 + 1. In our analysis we plot p^hr 1 ! 2 as a 
function of | In and vary the exponent S until one gets 
asymptotically a straight line. The best result is obtained 
for E = 0.41 and the corresponding plot is shown in 
Fig. | 
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p a (6p,h) = h 1 ' 2 |ln/i| E f(x,l), 



Similar to the lower dimensions we consider the scaling 
behavior of the order parameter as a function of the con- 
trol parameter for different external fields. In the inset of 
Fig. ^| we plot p a vs. p for different fields from h = 10 
up to h = 5 10~ 4 . Choosing h \- a IP\ lnA| s = 1 the scal- 
ing ansatz [Eq. (13f)] yields in leading order 
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FIG. 10: The scaling plot of the order parameter fluctu- 
ations Ap a for the four-dimensional model. The data are 
rescaled according to Eq.© using the mean-field value 
cr = 2. The inset displays the unsealed data, i.e., the fluctu- 
ations are plotted as a function of the density p for different 
values of the field h. 
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FIG. 11: The order parameter p a as a function of the field ft 
at the critical density p = p c for £) = 5. The dashed line 
corresponds to a power-law behavior according to Eq. (^) with 
P/o = 1/2. 

where the scaling argument x is given in leading order by 

x = Sph- 1 ' 2 \\nh\ b - s ' 2 (17) 

Varying the logarithmic correction exponents one gets a 
convincing data-collapse, which is shown in Fig.O, for 
£ = 0.49 and b - s/2 = -0.17. The first value is in good 
agreement with £ = 0.41 obtained from the scaling be- 
havior at the critical density [Eq. (|l5|)1 . Using the average 
value S = I + s/2 = 0.45 and b-s/2 = -0.17 we get the 
estimation £> = 6 + / = 0.28 which agrees with B = 0.24 
obtained from numerical simulations in zero-field. 

Furthermore we consider how the logarithmic correc- 
tions affect the scaling behavior of the fluctuations at the 
upper critical dimension. Similar to the order parame- 
ter [Eqs. ( p||l7| )] we assume for the leading order of the 
fluctuations the scaling behavior 

A Pa (S P ,h) = h-~ 1 ' /a \\nh\ v d(5ph- 1/a |lnfc| _, »,l) 

(18) 

with a — 2. It is known that the mean- field value of 
the fluctuation exponent is 7' = which corresponds to 
a finite jump of Ap a at the critical density 0, ||. In 
order to avoid that Ap a diverges for h — > at 8p = the 
logarithmic correction exponent T has to be set to zero 
too. Therefore we try to obtain a data collapse of the 
fluctuation data according to the ansatz 

Ap a (Sp,h) = JC^/i- 1 / 2 I ln/j,|-^, 1). (19) 

A good data collapse is observed for 77 = 0.39 and the 
corresponding plot is shown in Fig. [l^. 

V. ABOVE THE CRITICAL DIMENSION 

Above the critical dimension, i.e. D > 5, the scaling 
behavior of the CLG model is expected to obey again the 



FIG. 12: The scaling plot of the order parameter p a for the 
five-dimensional model, i.e., above the upper critical dimen- 
sion. The data are rescaled according to Eq. (|ll]) using the 
mean-field values /3 = 1 as well as a — 2. The inset displays 
the unsealed data, i.e., the order parameter p a is plotted as 
a function of the density p for different values of the field h. 
The dashed line corresponds to the zero-field behavior. 



scaling ansatz Eq. where the exponents are given by 
the mean-field values /? = 1, a = 2, and 7' = [|[ m, 
independently of the particular dimension. 

In Fig.[ll] we plot the field dependence of the order 
parameter at the critical density. For sufficiently small 
values of the field (h < 10~ 4 ) we observe the expected 
scaling behavior p a ~ h 1 / 2 . Thus the order parameter 
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FIG. 13: The scaling plot of the order parameter fluctu- 
ations Ap a for the five-dimensional model. The data are 
rescaled according to Eq. dli| ) using the mean-field values 
a — 2 and 7' = 0. The inset displays the unsealed data, 
i.e., the fluctuations are plotted as a function of the density 
p for different values of the field h. 
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TABLE I: The critical density p c and the critical exponents /3, a, V and 7 of the CLG model for various dimensions D. 
The values of the susceptibility exponent 7 are calculated via Eq. (IlOj). The symbol * denotes logarithmic corrections to the 
power-law behavior. The values of p c and (3 are obtained from [^|. 

D p c /3 <j Y 7 



2 0.34494 

3 0.21791 

4 0.15705 

5 0.12298 



0.637 ±0.009 
0.837 ±0.015 
1* 
1 



2.227 ±0.032 
2.075 ± 0.043 

2* 
2 



0.384 ± 0.023 
0.208 ± 0.042 

0* 




1.590 ± 0.033 
1.238 ± 0.046 
1* 
1 



p a .{8p 1 h) is determined for various fields in this regime 
(from h = 2 10 -6 up to 10~ 4 ). The curves are shown in 
the inset of Fig jl2| . Rescaling these curves according to 
the ansatz Eq. (gjwith 8 — 1 and a = 2 one gets a good 
data collapse (see Fig. |l2|) . 

Finally we consider the fluctuations of the five- 
dimensional CLG model. The inset of Fig. [l3| shows how 
the fluctuations are affected by the field h. With van- 
ishing field the curves become steeper until one gets a 
jump g I for h = 0. 

Considering the scaling behavior one has to taken into 
account that the jump corresponds to 7' = 0. Therefore, 
we plot the fluctuations Ap a as a function of 5ph~ 1 / 2 . 
The resulting data collapse is shown in Fig. |l3| and con- 
firms the assumed scaling behavior. 



VI. CONCLUSIONS 

We introduced a method which allows to apply an ex- 
ternal field in the CLG model. The external field obeys 
the particle conservation and is conjugated to the order 
parameter, i.e., it is realized as a spontaneous creation of 
active particles. We considered the order parameter as 
well as its fluctuations of the CLG model as a function 
of an external field in various dimensions (L> = 2, 3, 4, 5). 
Although the external field drives the system away from 
criticality the order parameter obeys certain scaling laws 
for sufficiently small values of the external field. These 
scaling laws are investigated and the corresponding expo- 
nents are determined numerically. The obtained values 
of the field exponent a are listed together with other crit- 
ical indices in Table |. At the upper critical dimension 
D c = 4 the usual scaling behavior has to be modified by 
additional logarithmic corrections. 
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VII. APPENDIX 

Recently, the scaling behavior of the well known Ising 
model was investigated at the upper critical dimension 
(d c = 4) and it was argued that the singular part of the 
free energy obeys the finite-size scaling ansatz || 

f L {t,h) = L- 4 /(ti 2 ln 1/6 L,hL 3 ln 1/4 L) (20) 

where t denotes the reduced temperature (t = T/T c — 1), 
h an applied magnetic field and L denotes the system 
size. Thus we assume that the free energy is a generalized 
homogeneous function 

f(t,h,L)= (21) 
XfitX- 1 / 2 ] In Al 1 / 6 , h\- 3 / 4 \ In AI 1 / 4 , LA 1 / 4 ) 

with A > 0. Of course for A = L~ 4 one recovers Eq. (|2"c|). 
This ansatz can be checked in the following way: the 
derivative of the free energy with respect to the applied 
field leads to the scaling equation of the magnetization 

m(t,h,L) = A 1/4 | In A| 1/4 (22) 
rh(t\- 1/2 \ In Al 1 / 6 , h\~ 3 / 4 \ In AI 1 / 4 , LA 1 / 4 ). 

Choosing <A _1 / 2 | In A| 4 / 6 = 1 one gets in leading order 
for the order parameter in the thermodynamic limit at 
zero field 

m(t, h = 0) -< 1/2 I lntl 1 / 3 . (23) 

The field dependence of the magnetization at the critical 
temperature (t = 0) 

m(t = 0,h) - h 1/3 |ln/i| 1/3 (24) 

is obtained by setting /iA _3 / 4 | In A| 4 / 4 = 1. Analogous 
one gets in leading order for the susceptibility and specific 
heat 

X (i,/i = 0)^t- 1 |lnt| 1 / 3 (25) 

and 

c(t,h = 0) - I ln^ 1 / 3 , (26) 

respectively. In this way the ansatz Eq. ( pi] ) leads directly 
to the Eqs. ( |2^j24] , |25| , p6| ) which were already derived in 
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the 1970s by Wegner and Riedel using renormalization 
group techniques JhJ. It is worth to mention that these 
are exact results within the renormalization group theory, 
i.e., neither the values of the mean-field exponents nor the 
values of logarithmic correction exponents are obtained 



from approximation schemes like e- or 1/n-expansions. 

In the case of the CLG model we choose a scal- 
ing ansatz for the leading order of the order parameter 
[Eq. (|l3|)] with corresponds to Eq. (p2|). 
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